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ABSTRACT 

We have entered the era of explosive transient astronomy, in which current and upcoming real-time 
surveys such as the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory {PTF) 
and Panoramic Survey Telescope and Rapid Response System {Pan-STARRS) will detect supernovae 
in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover 
supernovae from the earliest stars in the universe and reveal their masses. The observational signatures 
of these astrophysical transients are the key to unveiling their central engines, the environments in 
which they occur, and to what precision they will pinpoint cosmic acceleration and the nature of dark 
energy. We present a new method for modeling supernova light curves and spectra with the radiation 
hydrodynamics code RAGE coupled with detailed monochromatic opacities in the SPECTRUM code. 
We include a suite of tests that demonstrate how the improved physics and opacities are indispensable 
to modeling shock breakout and light curves when radiation and matter are tightly coupled. 
Subject headings: atomic data - methods: numerical - radiative transfer - supernovae: general 



1. INTRODUCTION 

We are now in the era of explosive transient astronomy, 
in which real-time all-s ky surveys w ith the Large Syn- 
optic Survey Teles cope (llvezic et aLlbOOS.) . the Palomar 
Transient Factory (|Law et al.ll2009[ ) and Pa noramic Sur- 
vey T elescope and Rapid Response System (jKaiser et al.l 
|2002[) are or soon will be detecting supernovae (SNe) in 
unprecedented numbers. Gamma-ray bursts now probe 
the universe out to z > 8 and upco ming telescopes 
such as the James Webb Space Telescope ([Gardner et al.l 
|2006[ ) could discover SNe from the earliest stars (Pop HI 
SNe) and reveal their masses. The observational signa- 
tures of these cosmic explosions are key to unveiling their 
central engines, the environments in which they occur, 
and to what precision they can be used as standard can- 
dles to pinpoint cosmic expansion rates and constrain the 
nature of dark energy. Light curves only provide indirect 
measures of the physics at work in transients but are by 
far the most abundant type of data. New advances in 
numerical codes can now produce SN light curves and 
spectra with greater accuracy and time resolution than 
in previous studies, mandating a new generation of sim- 
ulations to better understand these cosmic explosions. 

Numerous codes have been developed in the past 
three decades to simulate SNe and their light 
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curves and spectra, with varying degrees of physics. 
They include codes with radiation transport linked 
to simple hydro dynamics such as EDDIN GTON 
(lEastman fc Pintol J993) and others (Hocflic h et al.l 
|1993[ ), one- temperature radiation diffusion cod e s suc h 
as Kepler (jWeaver et al.l 119781: IWooslev et all I2002D . 
and radiation hydrodynamics (RHD) codes that 
evolve radiation and material t emperatures sepa- 
rately suc h as VISPHO T ( Ensman fe^Burrow^ [l992), 
TITAN d Cehmevr fc Mihalaa 11993) and STELLA 
(jBlinnikov et al. 1998). There are many codes which 
have been developed to post-process profiles from 
hydrodynamical co des with radiation tra nsport, in- 
cluding SEDONA (IKasen fc WooslevI 120091) . SYNOW, 
(Branch et al. '1985") and several Monte Carlo codes 
(Mazzah & Lucv 1993, Kronicr fc Sim 2009 and others). 
RHD simulations have also been post-processed to ob- 
tain non-LTE (local thermody namic equilibrium) spectra 
(jHauschildt fc Ensmanlll994f) . 

In this paper we present a new method developed at 
Los Alamos National Laboratory (LANL) in which we 
start from a full RHD simulation and post-process the 
resulting structure with detailed monochromatic opac- 
ities to calculate SN light curves and spectra with an 
unprecedented level of accuracy. This level of detail is 
necessary for times when radiation-matter coupling is 
important, from shock breakout until the start of the 
nebular phase when the SN ejecta become optically thin. 
Models with full RHD are required where radiation plays 
an important role in the hydrodynamics, such as where 
shock heating is important in determining the tempera- 
ture profile. The codes presented here do not incorporate 
non-LTE effects, which become more significant at late 
times in SNe. 

In the initial stages of supernovae, rebound from grav- 
itational collapse or thermonuclear burning drives a 
highly radiative shock into the upper layers of the star. 
The shock is not visible until it breaks through the sur- 
face of the star because opacity due to electron scattering 
in the intervening layers traps photons in the shock and 



they are advected outward by the flow. When it reaches 
the surface of the star, the shock abruptly accelerates 
in the steep density gradient there, becoming even hot- 
ter and releasing a sharp burst of photons into the sur- 
rounding medium. This results in an extremely bright 
breakout transient with peak bolometric luminosities > 
10^^ erg s~^, chiefly in the form of x-rays, hard UV and 
occasionally gamma rays. These radiative losses alter the 
hydrodynamics of the system and so cannot be ignored. 
The ambient medium into which the shock breaks 
out can strongly affect light curves and spectra in ways 
that may have been underestimated in past studies. 
The luminosity and hardness of the breakout transient 
are governed in part by how much the shock acceler- 
ates as it exits the star: the steeper the gradient, the 
brighter and harder the pulse. Radiation breakout can 
occur after shock breakout in dense circumstellar en- 
velopes because photons can be trapped by t he wind af- 
ter t h e shock exits the surfa c e of the star (lOfek et al.l 
\20Wt. IBalberg fc LoeE [2011 IChevalier fc IrwinI I2011L 
L.H. Frey et al., in preparation). Shock breakout 
has been the subject of nu merous analytical studies 
(iColgatd Il974t iMatzner fc McKcc 1999: Nakar & Sari 
120101 : iPiro et a l. 2010: Kat z et al] 120121) and numerical 
simulatio ns (lEnsman fc Burrows! 19921: iBlinnikov et al 



2000; To minaga et all 120091: iFrver et al 



201Qt iFrver et al.ll2010t iKasen et al.ll201ll) in the three 
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decades since it was first predicted (jColgati 119741 : 
iKlein fc Chevalien 119781 ). However, there have only 
recently been observations which are interpreted by 
some as shock breakout. These include direct ob- 
servations _o£X-ray pulses in GRB 060218/SN2006aj 
(ICampana 'erall^ [2006) and XRO 080109/SN2008D 
(jSoder bcrg et al. 2008), carly-timc peaks in UV light 
curves (Schawinski ot al. 2008; Gczari ct al. 2008), and 
indirect o bservations through IR echoes in the Gas A SN 
remnant (jDwek fc Arendtll2008D . 

The location and time at which radiation breakout oc- 
curs depend on the radius at which the star or its envi- 
ronment becomes optically thin and is a function of wave- 
length. Standard blackbody and light-crossing time ar- 
guments for peak breakout luminosities and pulse width 
assume a gray opacity and no attenuation beyond the r = 

I surface. In reality, the opacity can be a strong function 
of wavelength in SN ejecta and is poorly approximated 
by a gray opacity. Analyses of shock breakout obser- 
vations fi.e. iWaxman et al. 2007; Gheva lier & Fransson 
[200llGezari et al.ll2008l:[Soderberg et al.|[20 08^ have con- 
sidered a variety of progenitors and stellar environments 
but treat breakout as a wavelength-independent event 
that occurs at a single radius, temperature and time. De- 
pending on the wavelength band, this assumption can in- 
troduce significant errors to estimates of explosion prop- 
erties derived from observations. 

The breakout transient is generally followed by a rise 
and then gradual decline in bolometric luminosity that 
is powered by both the conversion of kinetic energy into 
thermal energy and the decay of ^^Ni. Shock heating 
from the supernova ejecta interacting with the surround- 
ing winds also contributes to the peak luminosity in Type 

II SNe, and to a lesser extent in Ib/c SNe. Shocks in 
Type la SN surroundings can also contribute to the light 
curves near peak luminosity. This luminosity powered 
by shocks and ^^Ni decay persists for several months in 



core-collapse and Type la explosions and for 2-3 years 
in pair-instability explosions. The fireball cools as it ex- 
pands and the peak of its spectrum evolves to longer 
wavelengths. At the same time, the envelope surround- 
ing the star that was ionized by the breakout pulse be- 
gins to recombine and blanket the spectrum with lines, 
particularly at shorter wavelengths. In the frame of the 
shock, the photosphere from which photons escape de- 
scends deeper into the ejecta as it expands. The heavy 
elements it exposes to the surrounding medium over time 
depend on how Rayleigh- Taylor instabilities and asym- 
metries mix the interior of the star as the shock prop- 
agates to its surface. When there is little mixing, lines 
due to heavier elements do not appear in the spectra un- 
til later times, but when mixing is extensive such lines 
can become manifest at early times. The order in which 
metal lines appear in SN spectra over time can be a pow- 
erful probe of mixing during the explosion, and of the 
central engine itself. 

Bolometric luminosities can surge at late times if 
the ejecta is dense and delays the escape of IR pho- 
tons or if the photosphere moves into a warm re- 
gion such as a ^^Ni layer. Dense winds can also im- 
pact spectra at late times by more thoroughly blan- 
keting the spectrum with lines, particularly at short 
wavelengths. The light curve can also brighten if the 
ejecta collides with a dense shell from a violent lumi- 
nous blue variable-type mass ejection (as in some Type 
Iln SNe). Such collisions can change even dim SNe 
into extremely bright e vents fe.g. [Smith fc McGravll2007l: 
van Marie et al.ll2010t [Moriva et al.ll2010l: iRoming et al.l 



20121 ). After breakout, a radiative precursor, which 



would not be seen in hydro-only codes, propagates in 
front of the shock and interacts with any surrounding 
material. If a shell was previously ejected, this precur- 
sor impacts the shell before the main shock and causes a 
bump in the light curve (D.J. Whalen et al., in prepara- 
tion). During the nebular phase, the influence of shock 
heating and radiative losses on the hydrodynamics are 
minimal and hydro-only codes can be used for accurate 
simulations. As departures from LTE become more sig- 
nificant during the nebular phase, the accuracy of the 
codes and opacities presented here will decrease. 

We describe the LANL RHD code RAGE and how it is 
used to model supernova explosions in Section [2] We re- 
view the LANL SPEGTRUM code, which post-processes 
RAGE profiles to compute time-dependent spectra and 
light curves in Section [31 and discuss how atomic opac- 
ities derived from the LANL OPLIB database are im- 
plemented in both codes in Section [31 In Section [51 we 
present the suite of tests used to verify SPEGTRUM. We 
analyze a series of tests in which more physical processes 
are progressively activated in RAGE in order to assess 
their effects on spectra and light curves in Section [SI We 
show why full RHD and accurate opacities are so im- 
portant to realistic supernova light curves and spectra in 
Section[7l and in Scction[8lwe summarize our conclusions 
and discuss future applications. 

2. RAGE 

RAGE (Radiation Adaptive Grid Eulerian) is a mul- 
tidimensional, multispecies Eulerian ada ptive mesh re- 
fincment RHD code developed at LANL (iGittings et al.l 
,2008. ) . It couples second-order conservative Godunov hy- 



drodynamics to either gray or multi-group flux-limited 
diffusion radiation transport. RAGE models radi- 
ation hydroflows on one-dimensional (ID) Cartesian, 
cylindrical and spherical polar coordinate meshes, two- 
dimensional (2D) Cartesian and cylindrical grids, and 
three-dimensional (3D) Cartesian boxes. The mass frac- 
tion of each species is evolved with its own continuity 
equation under the assumption that all species have the 
same velocity at a given mesh point. 

2.1. Microphysics 

Although RAGE has three-temperature (electron, ion, 
and radiation) capability, we use two-temperature (2T) 
physics in our simulations, in which radiation and mat- 
ter, while coupled, can be at different temperatures. 
This can occur during shock breakout, when matter 
is kn own to be out of t hermal equilibrium with pho- 
tons (jNakar fc Sarill2010[ ). This 2T physics is an im- 
portant improvement over the pure hydrodynamics and 
one-temperature radiation models used in previous SN 
light-curve calculations. However, deviations between 
the matter and radiation temperatures can imply de- 
partures from local thermodynamic equilibrium (LTE) 
and lead to inaccuracies in our opacities during break- 
out, as we discuss below. We do not explicitly transport 
gamma rays from the radioactive decay of ^^Ni in the 
ejecta. Instead, their energy is deposited locally in the 
gas according to 
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where TNi — 7.6x10^ s, tco — 9.6x10^ s and the mean 
energies released per atom for the decay of ^^Ni and^fCo 
are E n^ = 1.7 MeV and Eco = 3.67 MeV (Frvc r et al.l 
l2009f ). This approximation is valid at early times when 
the ejecta is dense and optically thick but can break down 
at later times when gamma-rays are no longer trapped. 
This occurs after about 50 days in Type la SNe and, de- 
pending on the degree of mixing, after 150 days in core- 
collapse SNe. The energy deposition described above 
has been tested by comparing it to non-thermal 7-ray 
transport, which p roduces the same result at early times 
(jFrver et al.ll2009| ). We apply the same constant speciflc 
heat Cy to each constituent material, so only ions con- 
tribute to the total pressure. We also use an ideal gas 
equation of state (EOS) because the densities commonly 
found in our astrophysical applications are off of existing 
EOS tables in RAGE by many orders of magnitude. 

2.2. Radiation Transport 

Both gray and multi-group flux-limited diffusion are 
implemented in RAGE with opacities derived from the 
LANL OPLIB database using the TOPS code, as ex- 
plained in Section S) Rosseland mean opacities are used 
to evaluate the diffusion coefficient in the transport equa- 
tions and Planck mean opacities are used to calculate the 
emission and absorption terms. All the tests in our study 
were done with gray flux-limited diffusion. OPLIB opac- 
ities are calculated under the assumption that the gas is 
in LTE and are thus functions of frequency, temperature 
and density only. At late times, when the SN becomes 
transparent, LTE may no longer be valid because indi- 
vidual fluid elements become too diffuse to remain in 



thermal and radiation equilibrium with their surround- 
ings. When this happens, the true opacity can deviate 
from those in our tables. The tests we describe in Section 
[6] show that LTE is a good approximation during shock 
breakout. We also include an advection term in the gray 
diffusion calculation. This advection term conserves en- 
ergy (consistently incorporates the work deflned by pres- 
sure and change in volume, P*dV) and includes mate- 
rial motion effects (assuming the radiation flows with 
the matter). Although this is not correct in the free- 
streaming limit, this error is likely to be smaller than the 
fact that we are using flux-limited diffusion to approxi- 
mate the transition from diffusion to free-streaming. 

2.3. Gravity 

Gravity can play an important role in core-collapse su- 
pernovae at very early times, causing ejecta to fall back 
onto the newly formed compact remnant. In most su- 
pernovae, the bulk of the fallback occurs in the first 1000 
s of the explosion (|Frver et al.l 1 19991 I2006al: I Zhang et all 
12008; Fryer 200^. In our calculations, we first evolve 
SNe until after fallback and nucle ar burning is complet e 
in codes s uch as Kepler, C ASTRO (lAlmgren et al.l|2010() . 
SNSPH (IFrver et al.l [20'06b) and others dPryerl [Tpa) . 
This is done in order to determine nucleosynthetic yields 
and how much material falls back onto the compact rem- 
nant so this material can be removed from the profiles 
prior to mapping them into RAGE. After fallback and 
burning is complete, gravity no longer plays a large role 
in the explosion and we neglect it in RAGE. As a conse- 
quence, the star can slightly expand into the surrounding 
medium in the time it takes the shock to reach the sur- 
face, but this motion is minor. If we regrid the explosion 
prior to breakout, the original proflle of the star is re- 
stored to the grid. 

2.4. Supernova Profiles 

Gas densities, velocities, speciflc internal energies (erg 
g~^), radiation energy densities (erg cm~^) and species 
mass fractions are required to initialize an explosion pro- 
flle in RAGE. Gas densities, mass fractions and veloc- 
ities are directedly imported from codes such as those 
listed above. However, rather than converting the gas 
energies in those codes into the speciflc internal energies 
required by RAGE, we calculate gas and radiation en- 
ergies from the temperature of the gas, assumed to be 
ideal and monatomic, which is also usually supplied by 
the codes: 
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Erad ~ a-T , 



(2) 



(3) 



where R is the ideal gas constant and a is the radiation 
constant. This avoids subtleties in how gas energies are 
deflned in various codes, some of which lump contribu- 
tions due to the ionization states of atoms and radiation 
in with the internal energy. In contrast to internal en- 
ergy, the temperature is unambiguously deflned in most 
codes. 

We map the proflle, which consists of the explosion, 
the surrounding star and the circumstellar environment. 



onto a ID spherical grid in RAGE. We choose the grid 
resolution so that the blast and its photosphere are well 
resolved, and we permit up to five levels of refinement 
in both the initial interpolation of the profile onto the 
grid and during the simulation. We set refiecting and 
outflow boundary conditions on fluid and radiation flows 
at the inner and outer boundaries of the mesh, respec- 
tively. When the calculation is begun, Courant times are 
initially small due to high temperatures, large velocities 
and small cell sizes. To reduce execution times and ac- 
commodate the expansion of the ejecta, we periodically 
regrid the profiles onto a larger mesh as the explosion 
grows. This significantly increases the time step on which 
the model evolves. We remap just the explosion itself, 
ignoring any medium beyond the shock or its radiation 
front, and then graft the original environment lying be- 
yond this radius back onto the shock or front on the new 
grid. We apply the same criteria in choosing a new grid 
as in the original setup. The entire calculation typically 
requires 20,000 CPU hours on LANL platforms. 

3. SPECTRUM 

We post-process RAGE hydrodynamical profiles with 
the LANL SPECTRUM code to construct light curves 
and spectra for our simulations. Stand-alone post- 
processing allows us to exploit the comprehensive LANL 
OPLIB opacity database to compute detailed absorption 
and emission features in our spectra without resorting 
to thousands of energy groups in the transport calcu- 
lations. These monochromatic opacities are calculated 
assuming LTE conditions for a range of temperatures, 
densities and materials; this database is described in dc- 
tafl in Section m SPECTRUM is general enough that it 
can accommodate hydrodynamic profiles from any code. 
The current version of SPECTRUM is parallelized and 
reads in pre-computed opacity tables, derived from the 
LANL OPLIB database, unlike the earlier version sum- 
marized in Fryer ct al. 2009. This previous version was 
serial, did not efficiently reference the frequency grid and 
required an additional external calculation to determine 
the opacities for each radial grid point. 

The calculations described below, with the monochro- 
matic OPLIB opacities in the spectrum calculations, in- 
clude the effects of line and continuum opacities, fiuo- 
rescence, Doppler shifting, time dilation and limb dark- 
ening. Other ways of calculating these effects, such as 
the Sobolev approximation, are not necessary since we 
include them explicitly. 

For ID RAGE simulations, densities, velocities, mass 
fractions and radiation temperatures from the finest level 
of refinement of the adaptive mesh grid are extracted and 
ordered by radius. These profiles typically contain more 
than 50,000 radial data points, but constraints on ma- 
chine memory and time prevent us from using all of them 
in a SPECTRUM calculation. A subset of the data from 
each RAGE profile must be carefully chosen to fully re- 
solve all radiating regions of the fiow from which photons 
can escape, such as the photosphere of an SN shock or 
the collision of a shock with a dense structure. The radial 
grid used in SPECTRUM for one set of SN simulations 
is described in more detail in Section [5.2.1l 

3.1. Luminosity of a Fluid Element 



The ID RAGE data, sampled as described above, is 
mapped onto a 2D grid of radial and angular bins in 
SPECTRUM as shown in Figure [U where the angular 
bins are uniform in /i = cos 9. SPECTRUM takes the 
fiuid at each mesh point to be in local thermodynamic 
equilibrium (LTE) with its surroundings and its radiation 
spectrum to be Planckian. In this approximation, each 
fluid element emits radiation at the same rate it absorbs 
it from its environment, so the Kirchhoff-Planck relation 
holds: 



f] — l^abst — Kabstiy, 



(4) 



where rj is the emissivity of the fluid parcel, Habs is its 
monochromatic absorption opacity at frequency v and 
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From Equations |4] and [5] we construct the luminosity of 
each fluid element into a given line of sight: 



U^^ = Kabs-m- 



2hv^ 
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(6) 



Here, r and /i deflne the radial and angular position of 
the fluid parcel, m, v and Tgas are its mass, velocity and 
temperature, v is the frequency at which it radiates in 
its own frame, and h^c and k are Planck's constant, the 
speed of light and Boltzmann's constant, respectively. 
Note that the spherical symmetry of the problem implies 
that the luminosity defined in Equation |6] is for the entire 
ring of material passing through the given mesh point in 
Figure [Hand centered on the fJ. — I axis. Consequently, 
m is the mass of the ring: 
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where pi is the ejecta density at r^ and ha is number of 
angular bins. 

Discretizing the spherically symmetric blast profile 
uniformly in fi guarantees that every fluid element at 
a given radius has the same volume and so the same lu- 
minosity. The v/c term in the denominator adjusts the 
luminosity to account for the relativistic time dilation as- 
sociated with the motion of the fluid and the v/c term in 
the numerator accounts for Doppler shifts in photon en- 
ergy due to motion relative to the line of sight. As shown 
in Figure [TJ if the ejecta is expanding photons from /i < 
are redshifted and from /i > are blueshifted. 

3.2. Escape Luminosity 

To calculate the luminosity from each fluid parcel that 
actually escapes into the interstellar medium we attenu- 
ate Lj/" ^ by e~'''»''j , where r^^^^^- is the integrated opti- 
cal depth from the mesh point out of the ejecta to pho- 
tons at frequency v. It is constructed from consecutive 
line segments as shown in Figure [2] 



'^ri.fj.j — / ^ ^ab 



iPi |rj+i 



(8) 




Fig. 1. — The SPECTRUM mesh. The radius of a given zone is defined at its lower face but /i 
uniformly partitioned in /i. 



cos 9 is cell-centered. The grid is 



where pi is the density along the segment (assumed to 
be that of the mesh point), and the sum is over all line 
segments extending from the source to the outer edge of 
the grid along the line of sight. As shown in Figure [H 
the Ti and /i^ of the fluid parcel uniquely determine the 
lengths of all segments from its position out to the edge 
of the grid because all the r are known and the ends of 
each segment are at the same height above the axis of 
symmetry. Each zone along the line of sight from the 
radiating cell intercepts photons at an energy different 
from that at the source that is determined by the total 
relative Doppler shift between them. Thus, the Kabs for 
each line segment is computed at the energy into which 
source photons arc Doppler shifted upon reaching that 
zone. 

Since Equation [5] is used to calculate the contribution 
to the escaping luminosity, the use of Kabs here is appro- 
priate as only the photons that are actually destroyed 
do not escape. However, in situations where there is a 
strongly scattering dominated atmosphere with little or 
no absorption, the thermalization depth of the photons 
will be significantly different from the depth of last inter- 
action. While these scattered photons are not destroyed, 
their properties may change due to a scattering event 
and their emergence from the ejecta will be undoubt- 
edly delayed in time. The approximation in Equation 
[8] will remain adequate for modeling the emergent spec- 
trum so long as two conditions are met: (1) the added 
delay from scattering processes near the thermalization 
surface {Tabs ~ 1) is short compared to the ejecta evo- 
lution timescale, and (2) the spectral information in the 
emerging photons is not strongly altered by the scatter- 
ing interactions. For typical conditions in low energy X- 
ray emission from supernova ejecta, Thompson electron 
scattering is sufficient and its coherent nature assures 
that the second criteria is met. In reality Thompson scat- 
tering limit is not perfectly coherent and attention must 
be paid to the number of scatters to ensure that the in- 
tegrated energy shift is not significant. The first criteria 
is one that must be checked for each new ejecta scenario 
and is an important caveat in the use of the SPECTRUM 
code. Scattering dominated atmospheres can occur in ex- 
tremely energetic explosions, such as the pair-instabilit y 
SNe of blue compact giant stars (jWhalen et al.|[2012aD . 
In these cases, the flux calculated in SPECTRUM with 
Ktot in Equation|S]instead of Kabs can be 10 orders of mag- 
nitude lower, with a peak wavelength 10 times longer, as 




Fig. 2. — The paths over which the T,-,,t are calculated. The 
ri and fij of a given fluid element set the lengths of all segments 
from its position out to the edge of the grid along the line of sight 
since the \r\ are simply the radii of adjacent mesh points and the 
ends of the segments all lie at the same height above the axis of 
symmetry. 

shown in Figure [3] 

Due to the high velocities seen at the shock front, even 
with well-resolved data there can be a significant veloc- 
ity gradient across just one zone. Consequently, photons 
from one end of the zone can be removed at a differ- 
ent energy at the other end that is determined by the 
difference in velocity between the two ends. Depending 
on the source frequency v, several energy bins may span 
the Doppler shift of the velocity gradient across the zone. 
We therefore frequency average the opacity of every zone 
by first computing its total opacity at each energy bin 
bracketed by the Doppler shift. We then sum them and 
divide by the number of bins across the shift to obtain 
the average opacity of the zone, effectively giving us sub- 
grid resolution for our opacities. The total luminosity 
escaping the explosion along the line of sight at source 
frequency v is then just the sum of the luminosities of 
each mesh point attenuated by the matter along the line 
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TABLE 1 
OPLIB Opacity Frequency Grid 
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Fig. 3. — Spectra from a 225 M© Pop III pair-instability SN at 
1.3 hours, calculated with either utot or K^bs in Equation [8] Note 
that luminosity is preferentially lost at short wavelengths. 



of sight 
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where na and ua are the number of radial and angu- 
lar bins, respectively. Equations [S] and IH] together with 
the symmetry of the grid ensure that limb darkening, 
Doppler shifts and time dilation due to relativistic ex- 
pansion of the ejecta are included in our spectrum calcu- 
lations. We also account for atomic emission lines with 
the Kabs term in Equation |6] and absorption lines with the 
g-'^'-i.fij term in Equation [51 Note that departures from 
LTE in the ejecta can lead to errors in the SPECTRUM 
opacities, and hence luminosities. Furthermore, devia- 
tions from LTE may also partially invalidate Equation |4] 
and cause additional errors in our spectra, especially at 
late times. 

The 14,900 energy bins for which there are opacities in 
the LANL OPLIB database, described in Section HTTl can 
change from point to point on the grid because the profile 
temperature varies with radius. These bins range from 
0.8 A - 2.0 X 10"^ A at 0.5 eV to 4.0 x 10"^ A - 100 A at 
10^ eV. To calculate the LIJ'^' one must choose a fixed 
set of 14,900 V for all the spectra. The bins we adopt are 
those for OPLIB opacities at T = 1 eV. These energies 
correspond to wavelengths of 0.415 to 9.1 x 10^ A, which 
capture all spectral features from shock breakout to the 
end of the light curve. We calculate Vj'^^ using the 
energy grid for the temperature in that cell and then 
map each L^, into the standard T — 1 eV energy grid. 

The SPECTRUM code can be extended to post- 
process data from 2D RHD simulations. The luminosity 
and the attenuation along a line of sight are calculated 
with the same methods described by Equations |6] and 
[51 The 2D discretization is done by taking parallel slices 
through the grid that align with the desired line of sight. 
Each slice can then be solved as separate ID spectra 
calculation and the total luminosity can be obtained by 
summing the spectra of individual slices. How we ex- 
tend SPECTRUM to general 2D data sets is described 
in detail in (W. Even et al., in preparation). 
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3.3. Light Curves 

Bolometric light curves can be constructed with SPEC- 
TRUM by summing over the full energy range of the 
spectra at each time step. Light curves for a wavelength 
band from a specific instrument can be created by con- 
volving the spectra and a filter response function. Each 
SPECTRUM calculation averages 200 CPU hours so the 
total CPU time required to compute a 250-point light 
curve is one to two times that of the RAGE run. How- 
ever, because many SPECTRUM jobs can run in simul- 
taneously and each one is completed quickly, it is usually 
possible to calculate a light curve per day on LANL plat- 
forms. 

4. OPACITIES 

The Rosseland and Planck mean opacities in RAGE 
and the monochromatic opacities in SPECTRUM are cal- 
culated with the LANL TOPS codel3 fro m cross s ections 
compiled in the OPLIB database (Ma gee et al.l [1995'). 
An advantage to using this database is that all of the 
opacities are calculated in a consistent manner, i.e. the 
results have not been culled from a variety of sources. 
A disadvantage is that only electric dipole transitions 
are considered in the bound-bound contributions, while 
other databases include forbidden lines. However, for 
LTE conditions. El-allowed transitions should be suf- 
ficient for capturing the relevant physics. We also note 
that the OPLIB database averages together some nearby 
spectral lines, but this level of accuracy is sufficient for 
the light-curve modeling considered here as well as for 
identifying general spectral features. 

We now describe these cross sections and opacities in 
detail and how mixed opacities are calculated for multi- 
species flows in RAGE and SPECTRUM. 



4.1. The LANL OPLIB Database 

TOPS can calculate external tables of absorption and 
scattering opacities k, for a set of density and temperature 
points that the user chooses. The k are derived from 
cross sections a in the OPLIB database that have been 
compiled for pure chemical elements that are assumed to 
be in LTE, according to 



i^u{p,T) 



TVo 
M 



<^:^{V,T), 



(10) 



where A^o is Avogadro's number, M is the atomic weight 
of the element and p is the density. The a are calculated 

^ http://aphysics2/www. t4.lanl.gov/cgi-bin/opacity/ tops. pi/ 
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Fig. 4. — The densities and temperatures for which hydrogen 
cross sections exist in the LANL OPLIB database. In general, the 
density range varies with both temperature and material. The red 
hatched lines indicate the original density range of the hydrogen 
cross sections and the black lines delineate our extension to this 
data set. 



from first principles based on quantum meclianics, and 
are tabulated on a grid of temperatures T, electron de- 
generacies r], and dimensionless reduced photon energies 
u given by 

hv , , 

This (T,77,u) grid is fixed for all elements. For a given ry 
and T, there are 14,900 u that range from 1.25 x 10~^ 
to 3 X IC*, which we summarize in Table [H 

TOPS can calculate gray, multi-group, or monochro- 
matic opacities for pure elements or mixtures of ele- 
ments and chemical compounds. Opacities in RAGE and 
SPECTRUM are referenced by p and T but, as men- 
tioned above, cross sections are tabulated by r\ and T in 
OPLIB. TOPS determines the p and T associated with 
an 77 and T according to the standard relation 



pkn.T) 



MNi{r],T) _ MN,{7j,T) 



No 
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(12) 



where p is the density of the pure element or mixture, 
Ni is its ion number density and Z is its average charge. 
The relation between N^ and rj is assumed to be that of 
a Fermi degenerate ideal gas: 



N, 



4 (27rTOefcT)3/2 



Fi/2{V), 



(13) 



where Fij2{r]) is the Fermi-Dirac integral of order 1/2, k 
is the Boltzmann constant, h is Planck's constant and ?Tie 
is the mass of an electron. When constructing an opac- 
ity table, TOPS will produce data only at temperatures 
that exist in the OPLIB database, i.e. neither interpola- 
tion nor extrapolation are performed for requested T that 
are not in the database. For this reason, the T points 
we adopt in our opacity tables are the same as those in 
OPLIB. On the other hand, if a requested density point 
falls between two corresponding p for which OPLIB a 



exist, TOPS will interpolate the cross section and com- 
pute an opacity for that point. If the density point falls 
above or below the corresponding range of p for which 
a are tabulated, TOPS returns the opacity of the clos- 
est corresponding on-table p and does not extrapolate. 
This procedure can lead to serious errors in opacity at 
very high or low densities that are beyond the OPLIB 
limits. The astrophysical simulations for which SPEC- 
TRUM was developed frequently have densities that are 
off-table on the low end of the original OPLIB grid. We 
therefore extended the database downward by 10 orders 
of magnitude in density for H, He, C, N and O, as shown 
in Figure 2] for H. The region in p and T for which opaci- 
ties of all the elements can be calculated is very similar to 
the red region shown for H; the extended regions for He, 
C, N and O are very similar to the black region show- 
ing our extension to OPLIB for H. The TOPS opacity 
tables are output to files that are read in by RAGE and 
SPECTRUM at execution time. 

We show how opacities for hydrogen vary with density 
and temperature in Figure [5] Absorption opacities at 
any given density generally fall as the temperature rises 
and the atoms become increasingly ionized (note that 
the total opacity falls to the Thomson scattering floor 
at sufficiently short wavelengths). Also, the absorption 
lines that are so prominent at lower temperatures fade at 
higher ones as the number of bound electrons decreases in 
the increasingly ionized gas. These trends demonstrate 
that the atomic physics included in the OPLIB database 
captures the effect of temperature and density on opacity. 

In certain applications it is advantageous to create ta- 
bles of opacities for a multi-material mixture with TOPS. 
To calculate mixed opacities, TOPS assumes that all 
constituents are homogeneously mixed in each cell (i.e. 
atomically mixed) and imposes the EOS condition of 
temperature-?7 (T-?^) equilibrium on the constituents of 
the mixture. By virtue of Equation 1131 imposing T-77 
equilibrium is equivalent to imposing T-A^g equilibrium, 
which is consistent with the physical picture of each con- 
stituent experiencing the same global free-electron distri- 
bution. Furthermore, since an EOS provides the relation 
between three physical quantities, imposing T-A^e equi- 
librium among the constituents is equivalent to imposing 
T-Pe equilibrium, i.e. the temperature and electron pres- 
sures associated with each constituent are equal to those 
of the mixture. Under these assumptions, TOPS obtains 
the opacity for a mixture at conditions p and T by sum- 
ming the cross sections for pure elements according to 



A^f 



Hu{p,T) = -Tj^fiCrlir],T), 



(14) 



where A'o is Avogadro's number, M is the atomic weight 
of the mix, fi is the atomic fraction of the i*^ element, 
and rj is chosen to satisfy Equation [12] TOPS calculates 
gray or multi-group mixed opacities by integrating these 
monochromatic opacities with the appropriate weighting 
function. 

4.2. RAGE Opacities 

We use TOPS to precompile Rosseland and Planck 
mean opacities for each pure element of the SN ejecta in 
a single external file for RAGE. However, since a mesh 
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Fig. 5. — Total opacity (upper row) and absorption opacity (lower row) for hydrogen as a function of temperature and wavelength. Left 
to right: T = 0.5 eV, 2.0 eV, and 100 eV. 



point on the RAGE grid might be populated by several 
elements, RAGE must construct a single opacity for the 
mixture. RAGE does not construct mixed opacities in 
the same way as TOPS because it imposes different EOS 
mix conditions on the constituents. Instead of taking 
them to be atomically mixed in each cell, RAGE assumes 
that they remain physically segregated (i.e. consistent 
with a heterogeneous or chunk mix), each with its own 
volume fraction and partial density (for additional de- 
tail s concerning the atomic and ch unk mix assumptions, 
see lSmith|[200l iPomranind fTOOl . Like TOPS, RAGE 
imposes T-P equilibrium on the components of the mix- 
ture, but considers both the ion and electron pressure. 
More specifically, in RAGE the ion pressure of each con- 
stituent is equal to the ion pressure of the mixture and 
the electron pressure of each component is the same as 
the electron pressure of the mixture. After obtaining the 
partial densities from the above procedure, the opacity 
of the mixture is obtained from 



s(p,T) = ^/fs,(p„r), 



(15) 



where k represents a mean (gray or group) opacity, /* is 
the mass fraction and pi the partial density of the i*-^ con- 
stituent, which is calculated from its volume fraction and 
mass fraction. The RAGE simulations presented in this 
paper use gray opacities, so R represents the Rosseland 
or Planck mean opacity. In regions where ion pressure 
is non-negligible and material is not atomically mixed. 
Equation [15] is expected to yield mixed opacities that are 



more accurate than those obtained under the assump- 
tion that el ectrons dominate the pressure, as in TOPS 
(jZeeb fc Po ntes 2006). However, because we are using 
an ideal gas EOS rather than the RAGE tabular EOS, 
we ignore electron pressure when calculating the pi. 

In principle, one should first construct the sum in 
Equation [15] at each frequency and then average over 
frequency to compute the mean mixed opacity rather 
than sum opacities that have already been frequency- 
averaged. In the Planck mean, the two approaches yield 
the same opacity because the mean is a simple arithmetic 
average for which summing and averaging are commuta- 
tive. The two operations are not commutative for Rosse- 
land opacities because the mean is harmonic, not arith- 
metic. Unlike TOPS, which frequency averages mixed 
monochromatic opacities to compute mean mixed opaci- 
ties, RAGE sums mean opacities, even though it is some- 
what less accurate, because computational constraints 
preclude a RAGE calculation for all 14,900 u. Finally, 
to construct the mixed opacity at a given mesh point at 
some p and T, RAGE performs a bi-linear interpolation 
of log R from Equation 1151 over log p and log T. 

4.3. SPECTRUM Opacities 

We use TOPS to create monochromatic absorption and 
total opacity tables for SPECTRUM, one for each ele- 
ment. We calculated tables for the first 30 elements at 
the 50 OPLIB temperatures from 5 x lO""* keV to 100 
keV and a grid of 81 densities from 10"^^ to 10^ g/cm^. 
Opacities for H, He, C, N, and O were calculated for an 



additional 20 densities from 10^^^ to 10^^^ g/cm'^, using 
data from the extended tables described earlier. Opaci- 
ties for all 14,900 u are tabulated for each p and T. The 
OPLIB u grid structure (shown in Table 1) simplifies re- 
trieval of the desired opacity from the tables because af- 
ter the u interval in which a given frequency v resides has 
been determined, a simple index that directly references 
that interval is easily computed. This approach results 
in spectrum calculation times that are much shorter than 
those done with bisection searches. 

To calculate the mixed opacity for a cell at a given p 
and T, SPECTRUM first computes the opacity Ki of each 
constituent of the mixture at its partial density pi. At 
the given frequency, log kJ, is extracted from the table at 
the two densities that bracket pi at the table temperature 
closest to T. SPECTRUM then performs a simple linear 
interpolation of the logarithm of these two opacities with 
respect to log pi while holding T constant to obtain the 
desired ni, i.e. as in TOPS, SPECTRUM does not in- 
terpolate kI in T. Like RAGE, the mixed opacity of the 
cell is then taken to be the sum of the opacity of each 
constituent weighted by its mass fraction, according to 

^,(p,r) = ^/f<(p„r), (16) 

i 

which is the monochromatic version of Equation 1151 
5. SPECTRUM TESTS 

We ran a series of tests to verify the SPECTRUM code 
and to determine the number and structure of radial bins 
and the number of angular bins required to converge to 
accurate spectra given memory and time constraints. We 
find that convergence is problem specific, and should be 
investigated for each new type of simulation. 

5.1. Blackbody Tests 

One basic test every spectra code should be able to 
pass is to reproduce an analytic blackbody in the case 
of a constant opacity. To do this, we calculated spectra 
for a static sphere with constant temperature, density 
and opacity, surrounded by a vacuum. To obtain accu- 
rate blackbody spectra, the opacity and density must be 
selected so that the r =20 surface is located a small dis- 
tance inside the outer surface. We initialized the sphere 
in SPECTRUM with a very fine radial resolution at the 
emitting surface. As we show in Figure [51 SPECTRUM 
accurately reproduces blackbody spectra for a range of 
temperatures. 

5.2. Convergence Tests 

As explained in Section [31 for ID problems SPEC- 
TRUM maps data onto a 2D grid with radial and an- 
gular bins. We test the convergence of our SPECTRUM 
calculation in both angle and radius with spectra from 
a fiducial 23 M© su pernova. The progenito r, which was 
evolved in TYCHO ([Young fc Arnettil2005[) . is a binary 
star that loses its H envelope at the base of the red giant 
branch, evolves as a Wolf-Rayet star, and then ends its 
life as a 6.4 Mp, CN O core that explodes as a Type Ic SN 
([Young et al.ll200"6l) . We first modeled the collapse of its 
core, the launch of the shock and all nuclear burning in 
a ID hydrocode ([Frveii[l"999( ) and then mapped its pro- 
file into RAGE together with a wind envelope, defined 
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Fig. 6.— Analytical (solid lines) and SPECTRUM (dotted lines) 
blackbody spectra at 1 eV and 100 eV. 

with a mass loss rate of 10~^ Mo/yr and a velocity of 10* 
cm/s. This star was originally investigated as a potential 
Gas A progenitor and is now being studied to understand 
the effects of pre-core collapse mass loss on shock break- 
out and SN light curves. The convergence tests presented 
here are for a single time after shock breakout, 400 s, but 
we have performed identical tests at a variety of times 
before and after breakout and find similar results. These 
radial and angular bin convergence tests were conducted 
in tandem over a 2D parameter space, 500-9000 radial 
bins and 2-400 angular bins, to find the best combina- 
tion of radial and angular resolution. 

5.2.1. Radial Convergence 

As noted in Section [3l the RAGE simulations used 
as input to SPECTRUM contain at least 50,000 radial 
data points but current memory constraints limit SPEC- 
TRUM input to less than 10,000 points, requiring the 
creation of a new radial grid. Densities, mass fractions, 
velocities, and temperatures from the finest mesh are ex- 
tracted from the RAGE profiles and sequentially ordered 
by radius. Three different binning schemes were con- 
sidered: linear, logarithmic, and a hybrid scheme which 
focuses resolution on the radiating region. A simple bin- 
ning scheme with uniform, linearly spaced, bins will com- 
pletely exclude the interior of the simulation and only re- 
solve the wind. Logarithmically spaced bins include the 
full range of radii but do not sufficiently resolve the ra- 
diating region around the shock. Figure [7| shows spectra 
for these two binning schemes. To obtain an accurate 
binning scheme, any region where rapid changes in den- 
sity, velocity or temperature occur and from which pho- 
tons can escape, such as around the shock or any shell 
or torus features in the wind, must be defined and suf- 
ficiently resolved in the new grid. The criteria used to 
define this new grid will change based on the details of 
each problem being studied. Generally the interior of the 
star, from which luminosity cannot escape, and the stel- 
lar wind, which is not finely structured, can be included 
at low resolution. A hybrid scheme which combines lin- 
ear and logarithmic binning was used for all subsequent 
spectra in this paper and is compared to simpler schemes 
in Figure [7] This scheme includes finely spaced linear 
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Fig. 7. — Output spectra for identical RAGE data from a Type 
Ic SN at 400 seconds after core-collapse, calculated with 160 angu- 
lar bins and 5200 radial bins organized by three different binning 
schemes: linear bins, log bins, and the hybrid scheme described in 
the text. 
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Fig. 8. — Radial convergence test, running SPECTRUM with 
the same RAGE data as [7] and 160 angular bins. This data is 
binned with varying radial resolution between the t = 20 surface 
and the radiation front. Note that the spectra calculated with 5000 
and 9000 radial bins are nearly identical. 



bins between the radiation front and the t—20 surface 
and more widely spaced log bins inside that radius and 
in the wind. The t—20 surface is calculated using the 
Thomson scattering opacity, and defines the minimum 
inner bound of the region from which luminosity can es- 
cape. The RAGE data inside each of these new radial 
bins is mass-averaged to ensure that very sharp features 
are included. 

To find the minimum number of bins that resolve the 
radiating region well, we performed a convergence test 
using between 500 and 9,000 radial bins in the region be- 
tween the T=20 surface and the radiation front, with the 
hybrid scheme described above. For each of these tests, 
the interior region and the wind are each resolved by an 
additional 100 logarithmically spaced bins. As shown in 
Figure m this SN Ic spectrum converges at about 5000 
radial bins. The spectra in this figure were all calculated 
with 160 angular bins. In addition to these SPECTRUM 
tests, we compared each binned parameter with the com- 
plete RAGE data set to confirm that the convergence 
shown here does correspond to a well resolved profile. 

5.2.2. Angular Convergence 

When the number of angular bins is increased in 
SPECTRUM, a finer grid is used to calculate luminosity 
along the line of sight for a single grid point (see Figure 
[TJ This results in a more accurate luminosity calcula- 
tion. A convergence study was performed with between 
2 and 400 angular bins and a fixed number of radial 
bins. Figure [9] shows that any number of angular bins 
greater than two gives a nearly identical continuum and 
that line intensities converge at 160 angular bins, with 
5200 radial bins. However, we note that higher resolu- 
tion might be required for detailed studies of individual 
lines. We again emphasize that convergence in general 
is problem-dependent because the optimum scheme for 
allocating grid points depends on the structure of the 
radiating flow. 

6. PHYSICS EFFECTS OF RADIATION HYDRODYNAMICS 
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Fig. 9. — Angular convergence test, running SPECTRUM with 
the same RAGE data as[7l 5200 radial bins and a varying number 
of angular bins. 

The numerical calculations of SN light curves per- 
formed to date have implemented varying degrees of 
physics, as noted in Section [2j To investigate the im- 
portance of radiation processes that have been neglected 
in past studies of SN explosions, we performed five ID 
simulations of the Type Ic SN described in Section 15.21 
in which additional physics was incrementally activated 
in RAGE. Because our goal is to examine how radia- 
tion physics affects SN profiles, spectra and light curves, 
we make no attempt to compare hydrodynamical algo- 
rithms. The five cases we consider are as follows: 

• 2T: Gray flux-limited diffusion (FLD) in which 
matter and radiation temperatures are evolved sep- 
arately in the simulation. This physics option is 
used in all the simulations presented in subsequent 
sections. 

• IT: Gray FLD in which matter and radiation are 
assumed to have the same temperature. 
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• 1ow-k: 2T gray FLD with k — 0.4 cm^/g (the sim- 
ple analytic scattering opacity for ionized H) for 
each constituent for all densities and temperatures. 

• high-K: 2T gray FLD with k = 10"* cm^/g for all 
densities and temperatures. Many astrophysics cal- 
culations have been done with codes that include a 
radiation pressure proportional to T"' in the EOS 
but do not perform radiation transport. Setting 
the opacity to a large constant value in RAGE can 
serve as a proxy for radiation pressure in the EOS 
because energy will be contained in the radiation 
field but not transported, except in extremely low 
density regions. This method allows gas tempera- 
tures to account for the presence of the radiation 
and prevents them from spiking to extremely high 
values in low density regions when internal energy 
is deposited in the cells (by shocks, for example). 

• hydro-only: no radiation transport is performed, 
just hydrodynamics, and there are no radiation en- 
ergies in the EOS. Because of its simplicity, the 
hydro-only option is useful for a wide variety of 
astrophysical problems in which radiation energies 
are negligible in comparison to internal energies. 
However, in problems such as SNe where regions 
of the domain are dominated by radiation, this 
method can lead to highly unphysical results. 

These models were evolved in RAGE to well past shock 
breakout, which occurs at ~ 0.7 s. 

The interior of the SN evolves differently in the five 
simulations, but this region is extremely optically thick 
and has little effect on the luminosity of the SN. The pho- 
tosphere of the shock, where gas abruptly becomes opti- 
cally thin, primarily determines the observational prop- 
erties of the SN. We plot gas velocities, temperatures 
and densities for the cell that has the maximum gas ve- 
locity on the grid, along with the position of this cell, as 
a function of time for all five tests in Figure [TUl 

We note that the IT and 2T simulation results are 
nearly identical throughout the duration of the IT simu- 
lation because radiation is still trapped. After radiation 
breakout, the explicit solver used by the IT simulation 
requires the use of extremely small time steps and it is 
impractical to run the simulation to completion. The 
hydro-only simulation has the fastest shock, with veloci- 
ties well above the speed of light. Shock temperatures in 
this model approach 1 GeV at breakout. These extreme 
velocities and temperatures arise because in a pure hy- 
drodynamic simulation kinetic energy from the shock can 
only be converted into internal energy when the shock 
breaks out into the low-density wind. In reality, this en- 
ergy should be partitioned between the internal energy 
of the gas Eint (oc T) and the radiation energy Erad 
(oc T^). In low density regions like the wind, most of 
the energy goes into radiation, and because the energy 
released by the shock is dominated by the T^ term and 
the gas and radiation temperatures remain coupled, the 
true gas temperature will be much lower. 

In the high-K test, velocities and temperatures in the 
cell with the highest gas velocity are consistently lower 
than in the hydro-only simulation. Although radiation 
is not efficiently transported in this model, it more prop- 
erly divides energy between gas and radiation when the 
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shock breaks out into the wind. The 2T and IT models 
are very similar to the high-K model prior to breakout. 
After breakout, the 2T and IT solutions still have similar 
velocities, but lower temperatures than the high-K case 
because the wind ahead of the shock is optically thin and 
radiation efficiently transports energy out of the shock 
and cools it. The sudden shift in the 2T model around 
50 seconds is an artifact due to the noisy velocity profile 
behind the shock. The maximum velocity is initially near 
the front of the shock but a higher velocity region forms 
slightly behind the shock, causing the maximum veloc- 
ity to suddenly shift to a deeper position. In the 1ow-k 
test, the cell with the highest velocity trails well behind 
those in the other simulations and remains at a higher 
density. In this run, the maximum velocity traces the 
expansion of the outer edge of the star because the low 
opacity prevents radiation from imparting momentum to 
the wind and accelerating it, as in the other models with 
radiation. 

This latter point can be more clearly seen in the tem- 
perature, density, and velocity profiles for the explosion 
in Figures 11111131 In the 1ow-k simulation, the radiation 
front is visible as the nearly uniform temperature region 
extending from the surface of the star into the wind at 
times greater than 0.1 s in Figure [TT] Comparing the po- 
sitions of the radiation front and the maximum velocity 
in Figure [T2I it is clear that the radiation has propagated 
well into the wind without depositing any momentum in 
it since its velocity is still the initial value of 10* cm/s. In 
the other four models, the position of the radiation front 
coincides with that of the shock that the radiation drives 
into the wind ahead of the expanding outer layers of the 
star. As shown in Figure I13[ this shock creates a more 
gradual transition in density between the wind and the 
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outer edge of the ejecta. The blowoff of the outer layers 
of the star changes the depth from which photons escape 
the flow. The behavior of the shock pr opagating into low 
densit y wi nd is examined in detail b y [Colgate fc Whitd 
(p6^ and lMatzner fc McKed jlpggl). 

We show spectra for the 2T, high-K, and 1ow-k models 
at 1 and 100 s in Figure UM They are calculated with 
SPECTRUM with the full set of monochromatic opac- 
ities described in Section 14.31 We exclude the hydro- 
only model because its 1 GeV temperatures yield spec- 
tra that are unphysical. The spectra follow the same 
trends as in Figure 1101 the emitting region is hottest 
in the 1ow-k simulations, followed by the high-K and 2T 
simulations. Bolometric luminosities for all three runs 
are shown in Table [2l They vary by up to seven orders 
of magnitude just because of differences in opacities in 
the radiation hydrodynamics (RHD) calculations. Their 
spectral profiles are also quite different, particularly at 
breakout. We list V-band and B-band magnitudes for 
all three models in Tables |3] and 21 Depending on the 
band and the time, any of the three simulations can be 
brightest, often at variance with their relative bolomet- 
ric luminosities. At early times when they are below 
the limits of detection, luminosities in both bands can 
vary by more than 100 magnitudes with model. At later 
times when they are above detection limits, luminosities 
in both bands can still vary by several magnitudes with 
model. These extreme variations in bolometric luminos- 
ity and visible magnitudes demonstrate the importance 
of accurate opacities and RHD to realistic SN light curves 
and spectra. Simulations with approximate opacities or 
no radiation transport can yield spectra with errors that 
are far greater than those in the observations. This test 
case demonstrates a subset of all the possible effects of 
not modeling the full physics of a SN. 

7. PHYSICS EFFECTS ON SPECTRA 

Analyses of shock breakout observations generally as- 
stune that the opacity is constant when deriving the ra- 
dius at which radiati on breaks free of the shock (i.e. 
iSoderberg et al.l 120081 ). In reality, radiation does not 
break free of the shock or wind at a single radius or time 
because the opacities of the outer layers of the star, ejecta 
or surrounding wind are all strongly dependent on wave- 
length, temperature and density. As shown in Figure 
[5l monochromatic opacities can vary by many orders of 
magnitude with wavelength, so a single radius and tem- 
perature cannot accurately describe this phenomenon. 
Figures [T51 \W[ [T71 and \TE\ all present data from the Type 
Ic SN progenitor described in Section [5^ 

To study radiation breakout, we calculated the r=l 
surface with monochromatic opacities for several wave- 
lengths as well as a constant ionized helium opacity, 
K=0.2. Figure [15] shows that for the wind of this SN 
opacity is high at short wavelengths, leading to a large 
r=l radius. After the shock breaks through the surface 
of the star and reaches the t = 1 surface for a given 
wavelength, the surface coincides with that of the shock 
thereafter and expands with it. As the ejecta expands it 
becomes spherically diluted and cools, causing the pho- 
tosphere to retreat into the ejecta at longer wavelengths 
at later times, as seen in the u-band after ^ 10^ s. 

Because radiation breaks free of the shock at a range 
of radii, luminosity estimates can vary widely with the 



choice of r = 1 surface, as we show in Figure [HI Here 
we plot bolometric luminosities derived from the Stefan- 
Boltzmann law, with the radius and temperature at the 
T = 1 surfaces shown in Figure [151 We contrast these 
semi-analytic light-curves with the bolometric luminos- 
ity derived from SPECTRUM calculations, which include 
luminosity emitted at multiple wavelengths and radii. It 
is clear that the bolometric light curves peak at very dif- 
ferent times depending on the choice of photosphere, and 
that the single-radius calculations dramatically overesti- 
mate or underestimate the luminosity. In the k—0.2 and 
u-band cases, the overestimates are partly because the 
Stefan-Boltzmann law fails to include flux attenuation by 
the envelope. The wind, primarily consisting of oxygen, 
has high opacities at X-ray wavelengths, which results 
in the large r =1 radii in Figure 1151 The temperature 
at these surfaces is low until the shock has propagated 
sufficiently far for the radiative precursor, and then the 
shock, to raise the temperature. This results in a bolo- 
metric luminosity from the X-ray r=l surface that is 
initially constant and then abruptly rises at ~ 2 x 10''' s, 
as the radiative precursor and then the shock heat the 
region outside the r =1 surface and increase the opacity. 
The late-time (after ^ 10** s) bumps in the k=0.2 and u- 
band light curves occur as the shock propagates outward 
and the r =1 surface for the given wavelength descends 
into hotter, more dense regions of the shock. The effect of 
radiation breakout occurring at different times and radii 
for various wavelengths is further illustrated in Figure 
1171 where we show light curves for six Swift ultraviolet 
and optical bands calculated with SPECTRUM. 

Some of the differences between the semi-analytic light 
curves and the SPECTRUM light curve in Figure [TBI are 
due to emission and absorption in the flow on either side 
of the shock, which cannot be taken into account by a 
simple Stefan-Boltzmann calculation. SPECTRUM cal- 
culates the luminosity that escapes the ejecta as a func- 
tion of both wavelength and radius, enabling us to study 
luminosity as function of both radius and time as shown 
in Figure [181 At 10, lO-* and 10^ s the double-peaked 
structure of the radiating shock is clearly visible. This 
feature is absent at 0.1 s because the shock is still in- 
side the star, and is smeared out at 10^ s because the 
ejecta has greatly expanded and become nearly trans- 
parent. Shortly after breakout, which occurs at '^ 0.7 s, 
luminosity primarily originates from the relatively nar- 
row shocked region itself but over time a radiative precur- 
sor propagates ahead of it, heats the material there, and 
causes it to radiate as well. As the ejecta expands and 
is spherically diluted, radiation can escape from deeper 
in its interior. This effect is visible in the profiles as a 
plateau that begins to extend from the rear surface of 
the radiating shock at 10'^ s. 

If all the photons broke free of the shock at once, the 
duration of the initial transient would be comparable to 
the light-crossing time of the star, since photons emitted 
from its poles and its equator would reach an observer 
at times that differ by the time it takes light to cross 
the star. However, radiation-matter coupling generally 
smears radiation breakout over much longer times than 
the light-crossing time of the star, as shown in Figure [TOl 
Here, we plot bolometric luminosities at shock breakout 
for a 225 Mq Pop III pair-instability supernova evolved 
in KEPLER and RAGE. Strong coupling between pho- 
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Fig. 15. — The radius of the t=1 surface vs. time, calculated at 
single wavelengths in the x-ray (10 A), m2 (2246 A) and u (3465 
A) bands, and with a constant helium opacity (k = 0.2). 

tons and matter causes the radiation front to blow off 
the outermost layers of the star with very high velocities, 
but this in turn impedes the escape of the photons into 
the surrounding medium. This delay is evident in the 
relative widths of the bolometric luminosities calculated 
by SPECTRUM and with the Stefan-Boltzmann law (we 
use KTh = 0.288, the opacity due to Thomson scattering 
in primordial H and He gas, to compute the radius of the 
r = 1 surface). The width of the breakout transient cal- 
culated analytically is roughly the light-crossing time of 
the star, given that its radius is a few x 10^'^ cm, but the 
breakout pulse calculated with RAGE and SPECTRUM 
is both broader and dimmer because photons remain cou- 
pled to the shock for longer times, as seen in Figures [T51 
andlHl 

8. CONCLUSION 

Current and future all-sky surveys will provide numer- 
ous light curves from supernovae and other transients, 
including many which do not fit a previously known 
type (for example, 3% of PTF SNe are unclassified as 
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helium opacity (k = 0.2), overlaid with the bolometric light curve 
from SPECTRUM. 

of March 20lf"l). However, light curves only provide 
indirect measurements of the physical properties of pro- 
genitor stars and their surrounding environments. Ana- 
lytical calculations performed to derive these properties 
assume a constant opacity and treat shock breakout as a 
single event. Approximations made in some hydro-only 
or RHD codes result in inaccurate opacities or incom- 
pletely coupled radiation and hydrodynamics, as demon- 
strated in Section |6l To correct some of the inaccura- 
cies inherent in these approaches, we have developed a 
new pipeline starting with RHD simulations in RAGE 
including accurate gray opacities. We then post-process 
these simulations with SPECTRUM, computing emitted 
and absorbed luminosity from accurate monochromatic 
opacities. This pipeline allows us to accurately model the 
luminosity from SNe as a function of radius, wavelength 
and time. 
Our RAGE and SPECTRUM tests show that shock 

-'-'' www.astro.caltech.edu/ptf/ 



16 






cn 
o 

E 




time (days) 

Fig. 17.— SPECTRUM light curves for six Swift UV/optical 
bands: v (5468 A), b (4392 A), u (3465 A), wl (2600 A), m2 (2246 
A), and w2 (1928 A). 



(D 



Fig 
with 



10^ 



1 




1 1 


1 1 

O.ls 

10s 

10^3 : 

lO^s 
lO^s 



10' 



10' 



10' 



radius (cm) 



. 18. — Luminosity as a function of emission radius calculated 
SPECTRUM, for five different times. 



breakout and radiation breakout in supernova explosions 
are inherently radiation hydrodynamical and wavelength 
dependent in nature, and cannot be accurately modeled 
by current analytical methods or numerical simulations 
that lack this physics. Doing so can lead to serious errors 
when deriving the properties of supernovae, and hence 
their progenitors and environments, from light curve ob- 
servations. We also find that inaccurate opacities can 
result in significant errors in light curves and spectra, 
even well after shock breakout. Another key to calculat- 
ing realistic SN light curves and spectra is ensuring that 
both simulations and postprocessing resolve the photo- 
sphere of the shock at all times, since this region of the 
flow is most important to its observational properties. 

The sensitivity of our tests to opacities and radiation 
physics and the detailed structure seen in spectra from 
monochromatic opacities suggest that even more accu- 
rate spectra can be obtained with multi-group radia- 
tion transport. However, because multi-group simula- 
tion times roughly scale as the number of groups, it is 
computationally prohibitive to run with a large number 




time (10 sec) 



Fig. 19. — Bolometric lightcurves, calculated analytically (dotted 
line) and with SPECTRUM (solid line), for the pair-instability 
supernova of a 225 Mq blue compact Pop III star. 



of groups. The results from a single-group RAGE sim- 
ulation post-processed with SPECTRUM can direct the 
choice of energy groups to cover the range of expected 
energies and have appropriate resolution near dominant 
spectral features. 

The pipeline we have developed with RAGE and 
SPECTRUM is now being used to model the observa- 
tional signatures of a broad range of SNe. Core-collapse 
supernovae in dense winds are being studied in which 
radiation breakout occurs in the wind, well after shock 
breakout from the surface of the star (L.H. Frey et al., in 
preparation). We are also now investigating supernovae 
with more complex environments, such as Type la explo- 
sions with circumstellar tori (W. Even et al., in prepara- 
tion) and Pop III core-collapse SNe in dense circumstellar 
shells (Type Iln SNe, whose brightness can rival those of 
much more massive pair-instability SNe, D.J. Whalen et 
al., in preparation). We have calculated light curves and 
spectra for both Pop III core-collapse and pair-instability 
SNe in order to determine out to what redshifts they can 
be detected hy JWST and f uture thirty-meter telescopes 
(jWhalen et al.ll2013l[20T2bl ia[). and are also modehng the 
explosions of supermassive Pop III stars in z ~ 15 pro- 
togalaxies, and pair-instability SNe in the local universe 
(C.C. Joggerst et al., in preparation). The flexibility of 
our approach will yield many insights into these astro- 
physical transients as well as future ones whose progeni- 
tors are yet to be characterized. 

We thank Gabriel Rockefeller of CCS-2 at LANL for 
the development of the code used to extract explosion 
profiles from RAGE AMR data. DJW was supported by 
the Bruce and Astrid Mc Williams Center for Cosmology 
at Carnegie Mellon University. Work at LANL was done 
under the auspices of the National Nuclear Security Ad- 
ministration of the U.S. Department of Energy at Los 



17 



Alamos National Laboratory under Contract No. DE- 
AC52-06NA25396. All RAGE and SPECTRUM calcu- 
lations were performed on Institutional Computing (IC) 



and Yellow network platforms at LANL (Conejo, Ma- 
pache, Lobo, Turing and Yellowrail). 



REFERENCES 



Almgrcn, A. S., Bcckner, V. E., Bell, J. B., Day, M. S., Howell, 
L. H., Joggerst, C. C, Lijewski, M. J., Nonaka, A., Singer, M., 
& Zingale, M. 2010, ApJ, 715, 1221 

Balberg, S. & Loeb, A. 2011, MNRAS, 414, 1715 

Blinnikov, S., Lundqvist, P., Bartunov, O., Noinoto, K., & 
Iwamoto, K. 2000, ApJ, 532, 1132 

Blinnikov, S. I., Eastman, R., Bartunov, O. S., Popolitov, V. A., 
& Woosley, S. E. 1998, ApJ, 496, 454 

Branch, D., Doggett, J. B., Nonioto, K., & Thielemann, F.-K. 
1985, ApJ, 294, 619 

Campana, S., Mangano, V., Blustin, A. J., Brown, P., Burrows, 
D. N., Chincarini, G., Cummings, J. R., Cusumano, G., Delia 
Valle, M., Malesani, D., Meszaros, P., Nousek, J. A., Page, M., 
Sakamoto, T., Waxman, E., Zhang, B., Dai, Z. G., Gehrels, N., 
Immler, S., Marshall, F. E., Mason, K. O., Moretti, A., 
O'Brien, P. T., Osborne, J. P., Page, K. L., Romano, P., 
Roming, P. W. A., Tagliaferri, G., Cominsky, L. R., Giommi, 
P., Godet, O., Kennea, J. A., Krimm, H., Angelini, L., 
Barthelmy, S. D., Boyd, P. T., Palmer, D. M., Wells, A. A., & 
White, N. E. 2006, Nature, 442, 1008 

Chevalier, R. A. & Fransson, C. 2008, ApJ, 683, L135 

Chevalier, R. A. & Irwin, C. M. 2011, ApJ, 729, L6+ 

Colgate, S. A. 1974, ApJ, 187, 333 

Colgate, S. A. & White, R. H. 1966, ApJ, 143, 626 

Dwek, E. & Arendt, R. G. 2008, ApJ, 685, 976 

Eastman, R. G. & Pinto, P. A. 1993, ApJ, 412, 731 

Ensman, L. & Burrows, A. 1992, ApJ, 393, 742 

Fryer, C. L. 1999, ApJ, 522, 413 

— . 2009, ApJ, 699, 409 

Fryer, C. L., Brown, P. J., Bufano, F., Dahl, J. A., Pontes, C. J., 
Frey, L. H., Holland, S. T., Hungerford, A. L., Immler, S., 
Mazzali, P., Milne, P. A., Scannapieco, E., Weinberg, N., & 
Young, P. A. 2009, ApJ, 707, 193 

Fryer, C. L., Colgate, S. A., & Pinto, P. A. 1999, ApJ, 511, 885 

Fryer, C. L., Herwig, F., Hungerford, A., & Timmes, F. X. 2006a, 
ApJ, 646, L131 

Fryer, C. L., Rockefeller, G., & Warren, M. S. 2006b, ApJ, 643, 
292 

Fryer, C. L., Ruiter, A. J., Belczynski, K., Brown, P. J., Bufano, 
F., Diehl, S., Pontes, C. J., Frey, L. H., Holland, S. T., 
Hungerford, A. L., Immler, S., Mazzali, P., Meakin, C., Milne, 
P. A., Raskin, C., & Timmes, F. X. 2010, ApJ, 725, 296 

Gardner, J. P., Mather, J. C, Clampin, M., Doyon, R., 

Greenhouse, M. A., Hammel, H. B., Hutchings, J. B., Jakobsen, 
P., Lilly, S. J., Long, K. S., Lunine, J. I., McCaughrean, M. J., 
Mountain, M., Nella, J., Rieke, G. H., Rieke, M. J., Rix, H.-W., 
Smith, E. P., Sonneborn, G., StiaveUi, M., Stockman, H. S., 
Windhorst, R. A., & Wright, G. S. 2006, Space Sci. Rev., 123, 
485 

Gehmeyr, M. k, Mihalas, D. 1994, Physica D Nonlinear 
Phenomena, 77, 320 

Gezari, S., Dessart, L., Basa, S., Martin, D. C, Ncill, J. D., 
Woosley, S. E., Hillier, D. J., Bazin, G., Forster, K., Friedman, 
P. G., Le Du, J., Mazure, A., Morrissey, P., Neff, S. G., 
Schiminovich, D., & Wyder, T. K. 2008, ApJ, 683, L131 

Gittings, M., Weaver, R., Clover, M., Bctlach, T., Byrne, N., 
Coker, R., Dendy, E., Hueckstaedt, R., New, K., Oakes, W. R., 
Ranta, D., & Stefan, R. 2008, Computational Science and 
Discovery, 1, 015005 

Hauschildt, P. H. & Ensman, L. M. 1994, ApJ, 424, 905 

Hoeflich, P., Mueller, E., & Khokhlov, A. 1993, A&A, 268, 570 



Ivezic, Z., Tyson, J. A., Acosta, E., AUsman, R., Anderson, S. F., 
Andrew, J., Angel, R., Axelrod, T., Barr, J. D., Becker, A. C, 
Becla, J., Beldica, C, Blandford, R. D., Bloom, J. S., Borne, 
K., Brandt, W. N., Brown, M. E., Bullock, J. S., Burke, D. L., 
Chandrasekharan, S., Chesley, S., Claver, C. F., Connolly, A., 
Cook, K. H., Cooray, A., Covey, K. R., Cribbs, C, Cutri, R., 
Daues, G., Delgado, F., Ferguson, H., Gawiser, E., Geary, 
J. C, Gee, P., Geha, M., Gibson, R. R., Gilmore, D. K., 
Gressler, W. J., Hogan, C, Huffer, M. E., Jacoby, S. H., Jain, 
B., Jernigan, J. G., Jones, R. L., Juric, M., Kahn, S. M., 
Kalirai, J. S., Kantor, J. P., Kessler, R., Kirkby, D., Knox, L., 
Krabbendam, V. L., Krughoff, S., Kulkarni, S., Lambert, R., 
Levine, D., Liang, M., Lim, K., Lupton, R. H., Marshall, P., 
Marshall, S., May, M., Miller, M., Mills, D. J., Monet, D. G., 
Neill, D. R., Nordby, M., O'Connor, P., Oliver, J., Olivier, 
S. S., Olsen, K., Owen, R. E., Peterson, J. R., Petry, C. E., 
Pierfederici, F., Pietrowicz, S., Pike, R., Pinto, P. A., Plante, 
R., Radeka, V., Rasmussen, A., Ridgway, S. T., Rosing, W., 
Saha, A., Schalk, T. L., Schindler, R. H., Schneider, D. P., 
Schumacher, G., Sebag, J., Seppala, L. G., Shipsey, I., Silvestri, 
N., Smith, J. A., Smith, R. C, Strauss, M. A., Stubbs, C. W., 
Sweeney, D., Szalay, A., Thaler, J. J., Vanden Berk, D., 
Walkowicz, L., Warner, M., Willman, B., Wittman, D., Wolflf, 
S. C, Wood-Vasey, W. M., Yoachim, P., Zhan, H., & for the 
LSST Collaboration. 2008, arXiv: 0805.2366 
Kaiser, N., Aussel, H., Burke, B. E., Boesgaard, H., Chambers, 
K., Chun, M. R., Heasley, J. N., Hodapp, K.-W., Hunt, B., 
Jedicke, R., Jewitt, D., Kudritzki, R., Luppino, G. A., 
Maberry, M., Magnier, E., Monet, D. G., Onaka, P. M., 
Pickles, A. J., Rhoads, P. H. H., Simon, T., Szalay, A., 
Szapudi, I., Tholen, D. J., Tonry, J. L., Waterson, M., & Wick, 
J. 2002, in Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, Vol. 4836, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, ed. 
J. A. Tyson & S. Wolff, 154-164 
Kasen, D. & Woosley, S. E. 2009, ApJ, 703, 2205 
Kasen, D., Woosley, S. E., & Heger, A. 2011, ApJ, 734, 102 
Katz, B., Sapir, N., & Waxman, E. 2012, ApJ, 747, 147 
Klein, R. I. & Chevalier, R. A. 1978, ApJ, 223, L109 
Kromer, M. & Sim, S. A. 2009, MNRAS, 398, 1809 
Law, N. M., Kulkarni, S. R., Dekany, R. G., Ofek, E. O., Quimby, 
R. M., Nugent, P. E., Surace, J., Grillmair, C. C, Bloom, J. S., 
Kasliwal, M. M., Bildsten, L., Brown, T., Cenko, S. B., Ciardi, 
D., Croner, E., Djorgovski, S. G., van Eyken, J., Filippenko, 
A. v.. Fox, D. B., Gal- Yam, A., Hale, D., Hamam, N., Helou, 
G., Henning, J., Howell, D. A., Jacobsen, J., Laher, R., 
Mattingly, S., McKenna, D., Pickles, A., Poznanski, D., 
Rahmer, G., Rau, A., Rosing, W., Shara, M., Smith, R., Starr, 
D., Sullivan, M., Velur, V., Walters, R., & Zolkower, J. 2009, 
PASP, 121, 1395 
Magee, N. H., Abdallah, Jr., J., Clark, R. E. H., Cohen, J. S., 
Collins, L. A., Csanak, G., Pontes, C. J., Gauger, A., Keady, 
J. J., Kilcrease, D. P., & Merts, A. L. 1995, in Astronomical 
Society of the Pacific Conference Series, Vol. 78, Astrophysical 
Applications of Powerful New Databases, ed. S. J. Adelman & 
W. L. Wiesc, 51 
Matzner, C. D. & McKee, C. F. 1999, ApJ, 510, 379 
MazzaU, P. A. & Lucy, L. B. 1993, A&A, 279, 447 
Moriya, T., Yoshida, N., Tominaga, N., Blinnikov, S. I., Maeda, 
K., Tanaka, M., & Nomoto, K. 2010, in American Institute of 
Physics Conference Series, Vol. 1294, American Institute of 
Physics Conference Series, ed. D. J. Whalen, V. Bromm, & 
N. Yoshida, 268-269 
Nakar, E. & Sari, R. 2010, ApJ, 725, 904 



18 



Ofck, E. O., Rabinak, I., Ncill, J. D., Arcavi, I., Cenko, S. B., 
Waxman, E., Kulkarni, S. R., Gal- Yam, A., Nugent, P. E., 
Bildsten, L., Bloom, J. S., Filippenko, A. V., Forster, K., 
Howell, D. A., Jacobsen, J., Kasliwal, M. M., Law, N., Martin, 
C, Poznanski, D., Quimby, R. M., Shen, K. J., Sullivan, M., 
Dekany, R., Rahmer, G., Hale, D., Smith, R., Zolkower, J., 
Vclur, v., Walters, R., Henning, J., Bui, K., & McKcnna, D. 
2010, ApJ, 724, 1396 

Piro, A. L., Chang, P., & Weinberg, N. N. 2010, ApJ, 708, 598 

Pomraning, G. C. 1991, Linear Kinetic Theory and Particle 
Transport in Stochastic Mixtures 

Roming, P. W. A., Pritchard, T. A., Prieto, J. L., Kochanek, 
C. S., Pryer, C. L., Davidson, K., Humphreys, R. M., Bayless, 
A. J., Beacom, J. F., Brown, P. J., Holland, S. T., Immler, S., 
Kuin, N. P. M., Dates, S. R., Pogge, R. W., Pojmanski, G., 
StoU, R., Shappee, B. J., Stanek, K. Z., & Szczygicl, D. M. 
2012, ApJ, 751, 92 

Schawinski, K., Justham, S., Wolf, C, Podsiadlowski, P., 

Sullivan, M., Steenbrugge, K. C., Bell, T., Roscr, H.-J., Walker, 
E. S., Astier, P., Balam, D., Balland, C., Carlberg, R., Conley, 
A., Fouchez, D., Guy, J., Hardin, D., Hook, L, Howell, D. A., 
Pain, R., Perrett, K., Pritchet, C., Regnault, N., & Yi, S. K. 
2008, Science, 321, 223 

Smith, C. C. 2003, Journal of Quantitative Spectroscopy and 
Radiative Transfer, 81, 451 

Smith, N. & McCray, R. 2007, ApJ, 671, L17 



Sodcrberg, A. M., Berger, E., Page, K. L., Schady, P., Parrent, J., 
Pooley, D., Wang, X.-Y., Ofck, E. O., Cucchiara, A., Rau, A., 
Waxman, E., Simon, J. D., Bock, D. C.-J., Milne, P. A., Page, 
M. J., Barentine, J. C., Barthelmy, S. D., Beardmore, A. P., 
Bietenholz, M. F., Brown, P., Burrows, A., Burrows, D. N., 
Byrngclson, G., Cenko, S. B., Chandra, P., Cummings, J. R., 
Fox, D. B., Gal- Yam, A., Gchrels, N., Immler, S., Kasliwal, M., 
Kong, A. K. H., Krimm, H. A., Kulkarni, S. R., Maccarone, 
T. J., Meszaros, P., Nakar, E., O'Brien, P. T., Overzicr, R. A., 
de Pasquale, M., Racusin, J., Rea, N., & York, D. G. 2008, 
Nature, 453, 469 

Tolstov, A. G. 2010, Astronomy Letters, 36, 109 

Tominaga, N., Blinnikov, S., Baklanov, P., Morokuma, T., 
Nomoto, K., & Suzuki, T. 2009, ApJ, 705, LIO 

van Marie, A. J., Smith, N., Owocki, S. P., & van Veelen, B. 
2010, MNRAS, 407, 2305 

Waxman, E., Meszaros, P., & Campana, S. 2007, ApJ, 667, 351 

Weaver, T. A., Zimmerman, G. B., & Woosley, S. E. 1978, ApJ, 
225, 1021 

Whalen, D. J., Frey, L., Even, W., Fryer, C., Heger, A., Holz, D., 
Hungerford, A., Lovekin, C, & StiavelU, M. 2012a, apj, 
submitted (arXiv:1211.4979) 

Whalen, D. J., Fryer, C. L., Holz, D. E., Hcgcr, A., Woosley, 
S. E., Stiavelli, M., Even, W., & Frey, L. H. 2013, ApJ, 762, L6 

Whalen, D. J., Joggerst, C., Fryer, C, StiavelU, M., Heger, A., & 
Holz, D. 2012b, arXiv:1209.5459 

Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Reviews of 
Modern Physics, 74, 1015 

Young, P. A. & Arnett, D. 2005, ApJ, 618, 908 

Young, P. A., Fryer, C. L., Hungerford, A., Arnett, D., 
Rockefeller, G., Timmes, F. X., Voit, B., Meakin, C, & 
Eriksen, K. A. 2006, ApJ, 640, 891 

Zecb, C. N. & Pontes, C. J. 2006, los Alamos National Laboratory 
Report No. LA-UR-06-0882 (unpublished) 

Zhang, W., Woosley, S. E., & Heger, A. 2008, ApJ, 679, 639 



